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We demonstrate that a Finite Bath Fluctuation Theorem of the Crooks type holds for systems 
that have been thermalized via weakly coupling it to a bath with energy independent finite specific 
heat. We show that this theorem reduces to the known canonical and microcanonical fluctuation 
theorems in the two respective limiting cases of infinite and vanishing specific heat of the bath. The 
result is elucidated by applying it to a 2D hard disk colliding elastically with few other hard disks 
in a rectangular box with perfectly reflecting walls. 
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I. INTRODUCTION 

During the last decade a number of fluctuation theo- 
rems have been reported in the literature, which have 
contributed a good deal to a better understanding of 
nonequilibrium thermodynamics [J, 0, S 0, S 0, HI- 
These can be roughly divided in two categories: steady 
state fluctuation theorems and transient fluctuation the- 
orems. The former apply to systems in nonequilibrium 
steady states and give information on the system fluctu- 
ations in the asymptotic regime of very large times (see 
O [HI for reviews on this topic) . The latter apply to 
systems that are temporarily driven out of equilibrium 
and give information about the fluctuations of work gen- 
erated by the driving forces. The most representative 
example of the latter kind is the Crooks fluctuation the- 
orem 0, ||| , that applies to systems that are initially in 
a canonical state. Although the canonical case is by far 
the most common case, one may need to study situa- 
tions where the system is initially distributed according 
to some other statistics, instead. For example the system 
might be initially a microcanonical state of well defined 
energy. In this latter case it has been shown that a mi- 
crocanonical fluctuation theorem of the type of Crooks 
also exists [H, [l3| . One naturally then wonders whether 
the same type of transient fluctuation theorem exists as 
well for yet other types of statistics. 

In this work we focus on the probability distribution 
function (pdf) that describes the statistics of a subsys- 
tem of a total classical ergodic system with fixed energy. 
For the case where the interaction between the subsystem 
and the rest of the total system (which we will refer to as 
the bath) is weak, and the specific heat of the bath is in- 
dependent of the energy (as for an ideal gas, or for a bath 
composed of hard spheres), the derivation of the pdf is 
a standard problem of statistical mechanics [HI, [HI . We 
make no assumptions regarding the size of the bath; in 
particular we do not assume that it is much larger or 
smaller than that of the system of interest as assumed 
in the canonical and microcanonical cases respectively. 



For this reason we refer to this type of bath as to a fi- 
nite heat bath, and to the statistics of the subsystem as 
to the finite bath statistics (see Eq. ([5]) below). For 
this statistics we show that a fluctuation theorem of the 
type of Crooks, i.e., a Finite Bath Fluctuation Theorem 
holds. This Finite Bath Fluctuation Theorem includes 
the Crooks canonical fluctuation theorem and the micro- 
canonical fluctuation theorem, as the two limiting cases 
in which the bath specific heat goes to infinity and zero, 
respectively. 

The present work is organized as follows. In Sec. [TTJ 
we review the derivation of finite heat bath statistics and 
recall some of its properties. In Sec. IHII we derive the cor- 
responding Finite Bath Fluctuation Theorem, and show 
that it reduces to microcanonical and canonical fluctua- 
tion theorems in the limits of vanishing and infinite baths, 
respectively. In Sec. |lV]we apply the theory to a specific 
example (i.e., a 2D hard disk elastically colliding with 
few other hard disks in a box) and test the validity of 
the Finite Bath Fluctuation Theorem, both analytically 
and numerically. Sec. [V] contains a discussion of the 
obtained results. The conclusions are drawn in Sec. ED 



II. FINITE BATH STATISTICS 

Let us consider a finite classical Hamiltonian system 
of total N to t particles and total energy E to t composed of 
two weakly interacting subsystems: the "system of inter- 
est" (or simply the system) and the "bath". Assuming 
that the total system is ergodic, the probability density 
function (pdf) of the system is given in terms of density 
of states, ftB(E), of the bath and density of states of the 
total system, H, to t(E), as (l6| : 



p(z, A) 



n B {E tot -H(z,X)) 

Qtot(Etot) 



(1) 
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where z = (pi, p s , qi, q s ) stands for the 2s dimen- 
sional phase space point of the system. Here we assume 
that the (sub)-system Hamiltonian H(z, A) may depend 
on some externally controllable parameter A (this could 
be for instance the volume of a vessel that contains the 
system, or an applied magnetic or electric field). For ex- 
ample, in the case of a bath composed of n hard spheres 
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in 3 dimensions, it is Qb(Er) oc E^ 1 ^ 2 1 (see Appendix 
IS]) , and one finds from (JT|) [14| : 



p(z, A) 



[Etot ~ H(z, A)] 



3n/2-l 



/(fa^-HCz.A)]^ 2 - 1 



(2) 



which is a known result of classical statistical mechan- 
ics [l7|- The symbol [a;]+ is defined as [x]+ ■= x6(x), 
with 6*(x) denoting Hcavisidc step function. Note that, in 
this case, the specific heat of the bath C(Eb), is energy- 
independent and equal to 3n/2 flSll . More generally one 
has the following theorem [Tg, 13 : 

Theorem 1 The system pdf is given by 
[E tot - H(z, A)] 



p{z, A) 



c-i 
+ 



/ dz[E tot -H(z,X)} 



c-i 



(3) 



z/ and onZy «/ t/ie specific heat of the bath C is energy- 
independent. 

Here C is the microcanonical specific heat of the bath, 



C(E E 



dEj 



-Tb(Eb) 



(4) 



where Eb is the energy, and Tb{Eb) '■= 
®b(Eb)/Qb{Eb) is the microcanonical tempera- 
ture expressed in terms of the phase space volume 
^b(Eb) of the bath, with energy below Eb, and the 
bath density of states VL B {E B ) = d$B(EB)/dE B - In the 
following of this work we restrict ourselves to the case of 
energy independent, positive specific heat of the bath, 
abbreviated as C(E B ) ■= C > 0. 

Note that the pdfs in Eq. ([3]) are parametrized via 
the total system energy E to t ■ It is however convenient to 
parametrize the pdfs via a property that pertains to the 
subsystem only, e.g., its average energy U. This is accom- 
plished by writing E to t = U + CT (here CT represents 
the average energy of the bath) , substituting this expres- 
sion in ([3]), and imposing that U = J dzH(z; X)p(z, A). 
This leads to solving the following equation for T, given 
the average energy U and A 



/ dzE(z- A) [1 - (H(z; A) - U)/(CT)} 



c-i 



Jdz [l-(H(z;X)-U)/(CT)} 



c-i 



U (5) 



We shall denote the value of T that satisfies Eq. §5§ 
for given U and A as T(U, A) (in Appendix [B] we prove 
that a solution T(U, A) always exists). With this function 
at hand we can parametrize the pdfs in Eq. ([3]) via the 
subsystem average energy U and recast them in the form: 



c-i 



, TT u [l-{H{z-\)-U)/(CT(U,X))] + 

PC{Z] U > A) = Nc(U^\) 

where Nc{U, A) is the normalization: 

Nc(U, A) = fdz [1 - (H(z: A) - U)/(CT(U, A))] 



(6) 



c-i 



(7) 



As discussed in Appendix [B] it is not always possible 
to invert T(U,X). For sake of simplicity, in this work 
we shall assume that T(U, A) is invertible with respect to 
the argument U. This means that we could also choose 
T as an independent parameter and express U as a func- 
tion of T and A. Thus we are free to choose between two 
possible parameterizations: a microcanonical- like param- 
eterization (or [/-parameterization), and a canonical-like 
parameterization (or T-parameterization) (2l| . 

We shall refer to the numerator in Eq ([6]) as to a "gen- 
eralized Boltzmann factor" . It is important to stress 
that a factor of the type [1 - (H(z; A) - [/)/(CT)]^ _1 
is a generalized Boltzmann factor only if T = T(U, A), in 
agreement with Eq. (|5]). 



Remark 

By expressing the specific heat C as C = 1/(1 — q) > 
one recognizes that the pdf in Eq. is the Tsallis 
escort pdf of index q with q < 1 [22|]. Note that these 
do not exhibit heavy tails but rather have a faster than 
exponential decay with a finite cutoff occurring at the 
energy U + CT = Etot- The physical meaning of this 
cutoff energy is that the system's energy cannot be larger 
than the total energy. 



Properties 

Equipartition 

The following equipartition theorem holds for the Fi- 
nite Bath Statistics in © [H|: 



(8) 



where (•) denotes average over pc in Eq. ([6]), pi is one of 
the momenta and repeated indices are not summed. Eq. 
© says that T(U, A) can be interpreted as the tempera- 
ture of the system. 



Heat Theorem 

The Finite Bath Statistics provides a mechanical model 
of thermodynamics [23| . meaning that the temperature 
T, the external parameter A, its conjugated generalized 
force f\, and the average energy U are related in such a 
way as to satisfy the heat theorem [24j : 



dU + f x dX 



exact differential 



where f\ is defined in the usual way as: 

f - ( dH ) 
fx -(^) 



(9) 



(10) 
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This property is an important one because it allows to 
determine the thermodynamic entropy associated with 
the Finite Bath Statistics by finding the integral of the 
exact differential. This is given by [24| : 



Sc(U,X)=\nN c (U,X) 
Interpolation 



(11) 



The pdfs in Eq. ^ interpolate between canonical and 
microcanonical ensembles. Using the limits of infinite 
and null specific heat C, i.e., 



limc^oo 


H 


X ' 

' c. 


C-l 

+ 


= lim 

C— oo 


H 


X ' 

' C- 


° = e x ; (12) 

+ 


lime— o 


H 


X ' 

' c. 


c 
+ 


= 6(x); 






(13) 


lim c _ 


H 


X ' 

' c. 


C-l 

+ 


= S(x) 






(14) 



one recovers the canonical and microcanonical pdfs [25j | : 

e -H(z;A)/T 



lim p c (z;T,A) 

C— *oo 



Z(T, A) 



(15) 



respectively [2(J. The microcanonical normalization 
r2(t7, A) is the system density of states. Likewise one has, 
for the normalization, the following limits [20| : 

lim N C (T,X) = [ d Z e-( H ^- u ^ T 

C— oo J 

= e u/T Z(T,X) (17) 
dz = $([/, A) (18) 



lim N C (U,\) 
c— o 



H{z;\)<U 



The quantity $([/, A) is the volume of system phase space 
with energy below U . The density of states is related 
to $ via a partial derivative f2 = d^/dU. By taking 
the logarithm one recovers canonical and microcanonical 
entropies; i.e., 

lim S C (T,X) = ^+lnZ(T,A) (19) 



Yrni^S c {U,\) = ln$([/,A) 



III. THE FLUCTUATION THEOREM 



(20) 



Consider an ensemble of systems distributed accord- 
ing to Eq ©. Assume the system being decoupled from 
its bath and that it is acted upon by an external force 
that changes the external parameter A according to some 
prescribed protocol \(t) executed between times to and 
t j . The probability density that the external force does 



a certain work W on the system in that interval of time 
reads: 



J dz Q 5(H f (z f )-H (z )-W) 



H a (z a ) - U 
CT (U) 



c-i 



where Zf = z(tf ,toi z o) is the solution of Hamilton's 
equation with initial condition zq. For simplicity of nota- 
tion we drop the variable A in all quantities that depend 
on it, and replace it with a subscript or /, depending on 
whether the quantity is taken at values of A equal to X(to) 
or X(t f ), e.g., H (z) = H(z,X(t )), T (U) = T(U,X(t )). 
By making the change of variables from Zq — > Zf with a 
unitary Jacobian, one obtains 

Nc,o(U)p?; U to (W) = I dz f 6(H (z ) - H f (z f ) + W) 



H f (z f )-(U + W) 



n C-l 



CT (U) 



(22) 



where now Zo = z(to, tf, zj), is the solution of Hamilton's 
equation with as initial condition and time running 
backward. Note that the second term in the integrand 
is not a generalized Boltzmann factor because in general 
it does not satisfy Eq. ©. However for any ST one can 
rewrite the previous equation as: 



n c.u m/ x (T {U) + 5T 



c-i 



J dz f 6(H (z a )-H f (z f ) + W) 
H f {z f ) - (U + W -C5T) lC 



1 - 



C(T (U)+6T) 



(23) 



We now choose 5T as the solution of the following integral 
equation: 



/ dzH f {z)B(z,U,W,ST) 



U + W- C8T 



J dzB(z, U, W, ST) 
where, for convenience, we introduce the notation 



(24) 



B(z,U, W 7 ST) 



1 



Hf(z) — (U + W — C8T) 



c-i 



(25) 



C(T (U)+ST) 
or, equivalcntly as a solution of: 

T (U) + ST = T f (U + W -C5T) (26) 



4 



Then, wc find: 

N c ,o(u)p?; u to ( w ) 



T f {U + W 



To(U) 



C6T) \ 



c-i 



dz f S(H (z )-H f (z f ) + W) 

H f (z f ) - (U + W -C5T) 
T f (U + W -C8T) 



c-i 
(27) 

where the second term of the integrand is the Boltzmann 
factor of the pdf pc( z 7 U + W — C5T, X(t f ) ) . The integral 
is the product of N C j(U + W -CST) and the probability 
p? '%+ W ~ C6T (-W) that the force performs the work -W 
when the protocol is run backward and the system is 
initially in the state p c (z: U + W - C5T, X(t f )). 

Therefore the following fluctuation theorem is ob- 
tained: 

To, 



c,u f 



N c ,o(U) ' 



(28) 



where, 



U f := U + W-CST (29) 

Tf := T f {Uf) (30) 

Using Eq. (TlTj) , Eq. (J25J) can be rewritten in terms of 
entropy as: 

c-i 



p?; u t0 (w) 
pf£(-w) 



exp[AS%°(U,W)] (31) 



where AS&°(U,W) = S CJ (U f ) - S c .o(U). 

The Finite Bath Fluctuation Theorem of Eq. J3II al- 
lows to calculate the ratios of probability of work done 
on the system when it is driven arbitrarily away from 
equilibrium during the action of the forward and back- 
ward protocol, in terms of equilibrium properties such as 
entropy and temperature. 



Alternatively one can take the limit C — > 
of Eq. If2"8"|) directly and obtain the expression 
T (U)$ f (U + W)/(T f (U + W)$o(U)), which reduces to 
the previous one by virtue of the microcanonical equipar- 
tition theorem. 



Limit of canonical ensemble 

Likewise, using the T parameterization, it can be seen 
that, in the limit C — > oo Eq. ((22)) becomes 

Z (T)p?;l(W) = e w / T 

x J dz f S(H (z ) - H f (z f ) + W)e- Hf{zf)/T 

(34) 

One thus obtains the fluctuation theorem for the canon- 
ical ensemble of Crooks 0, Q : 



pT;J (w) _ Zf (T) 



pTJA-w) z (t) 



W/T 



(35) 



IV. EXAMPLE: A 2D GAS OF HARD DISKS 

In this section we illustrate the Finite Bath Fluctua- 
tion Theorem by applying it to a system composed of 
n + 1 clastically colliding hard disks in a 2-dimensional 
box with perfectly reflecting walls. One disk will be our 
system of interest, whereas the remaining n ones will 
form the bath. We assume that the disks do not have ro- 
tational degrees of freedom. As shown in the Appendix 
IA1 the specific heat is given in this case by C = dn/2 
where d is the number of translational degrees of free- 
dom of each disk. In this case d = 2, hence C — n. Note 
the fact that C does not depend on energy 



Recovering known special cases 

Limit of microcanonical ensemble 

In the limit C — > Eq. ([22]) becomes (using the for- 
mula 6 (ax) = a~ 1 (5(a;), and Eqs. dHJ) and (fig)) ) 

t> (U)p?; u to (W) = T (U) 



x J dz f 5(H (z ) - H f (z f ) + W)S(H f (z f ) - (U + W)) 

(32) 

Using the microcanonical equipartition theorem [l6j 
T(U, A) = $(U, X)/fl(U, A), one recovers the microcanon- 
ical fluctuation theorem 113, 03 : 



p°t} u t0 (w) _ n f (u + w) 



o,u+w ( 

PtQ.tf \ 



-W) 



n (u) 



(33) 



The probability density function 

The energy of the system of interest is simply its kinetic 
energy; i.e., 



H(p x ,p y ;M) 



pI+pI 

2M 



(36) 



which fluctuates permanently due to the collisions with 
the bath's particles. According to Eq. ©, the probabil- 
ity that the disk has a given momentum (p x ,Py) is given 
by 



pc(px,p y ;U,M) 

= Nc\U,M) 



( P 1+pI)/(2M)-U 
CT(U, M) 



-\ c-i 



(37) 
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We consider the mass of the disk M as an external param- 
eter that can be changed at will in the course of time ac- 
cording to pre-spccificd protocols. The function T(U, M) 
has to be computed via Eq. ((SJ). In general, the solution 
of Eq. ([5]) with a purely kinetic Hamiltonian with s trans- 
lational degrees of freedom gives the usual equipartition 
of energy [22|: T(U,M) = 2U/s. In the specific case of 
Eq. HMD s = 2, hence 



T(U, M) = U , 



(38) 



and 



pc(px,p y ;U,M) 
= Nc 1 (U,M) 



1 - 



Vl +P 2 y 

2M 



-1 c-i 



- u /(CU) 



Using Eq. ©, with d39j) gives 

N C (U, M) = 2ttA[\ + C~ llC 



MU 



(39) 



(40) 



where A is the reduced volume (i.e., area in this 2- 
dimensional case) of the box (see the Appendix [A] for 
the definition of reduced volume). From Eq. (|3"9")l . one 
obtains the pdf of energy E of the disk: 



p(E;U) = U~ l [l + C 



-li-C 



1 - 



(E-U ) 
CU 



n C-l 



(41) 



Interestingly, the energy pdf does not depend on the mass 
M. In Fig. [T]we compare Eq. (|4"Tj) . with the result of 
various numerical simulations with C = 1,2,3,4. Note 
that for C = 1 the distribution is flat, for C = 2 it is lin- 
ear, for C — 3 it is quadratic etc... In view of theorem 1, 
the impressive agreement between theory and numerics 
corroborates the validity of the assumed ergodic hypoth- 
esis for this model system. Similar simulations have been 
reported in [2o] | for a 1-dimensional harmonic oscillator 
coupled to a bath of n 1-dimensional quartic oscillators. 
In that case the density of states of the bath is propor- 
tional to £;(3«-2)/4 j and 

accordingly the specific heat, 
C = (3n + 2)/4, is energy independent. 



Analytical test of the Finite Bath Fluctuation 
Theorem 



Consider a protocol M(t) that changes the mass of the 
disk from the value M = M(t ) to M f = M(tf). Ac- 
cording to the general assumption of our derivation, the 
system is decoupled from the bath during the action of 
the protocol. We are interested in checking the validity 
of Eq. (|28p . To this end we need to compute the for- 
ward pdf of work, Ptj,t (W), the backward pdf of work 

Pt^ti{—W), and the starting average energy of the back- 
ward protocol Uf, given the starting average energy of 
the forward protocol U . Solving Eq. (|2T))) with Eq. 




1.5 2.0 
E (kJ/mol) 

FIG. 1: (Color online) Energy probability density function 
(pdf) for a 2D hard disk of radius r = 1 nm and mass M = 2 
amu, in a bath composed of 1 (•), 2 (■), 3 (♦), 4(a) other 
identical disks. The dots represent histograms of properly 
normalized relative frequencies from numerical simulations. 
All simulations were carried out for the same total energy 
Etot — 3.3469 kJ/mol, which corresponds to measured av- 
erage energies of the disk of interest Ui = 1.67166 kJ/mol, 
U 2 = 1.11644 kJ/mol, U 3 = 0.835784 kJ/mol, U 4 = 0.671521 
kJ/mol. The solid lines represent the pdf predicted by the the- 
ory (Eq. ijHJ) for the measured average energies Ui,i = 1..4. 



(note that Eq. (|38p does not depend on the value of M, 
hence T f (U) = T (U) = U) wc arrive at: 



ST = W/(1 + C) 
hence from Eq. ([2^]) we obtain 

U f = T f = U + W/(1 + C). 



(42) 



(43) 



Using Eq. (|40[) with (|4"3"|) we obtain the normalizations 
of the equilibrium pdfs with average energy and external 
parameters (U,Mq) and (Uf,Mf) respectively: 



N c ,o(U) = 2ttA[1 + C- 1 ] l 'M U 



-1]Cj 
-llC; 



N C j{U f ) = 2TrA[l + C- L ] u M f [U + 



W 
l + C 



Using (I43I44|45|) we find: 

(TfY' 1 N CJ (Uf) 
\T J N c ,o(U) 

From Eq. (f2"Tj) we have: 
P?; U t (W) :=Ncl(U)A 

/ 2 

dp x dp y 8 



Ut\ C M 



f 
Mo 



(44) 
(45) 



(46) 



Pi 



Px+Pv 



2M f 



2M a 



- W 



Px+Py 

2M 



- U /(CU) 



c-i 



(47) 
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-1.4 -1.2 -1.0 -0.8 -0.6 -0.4 

W (kJ/mol) 



0.2 0.0 



FIG. 2: (Color online) Comparison between the numerical 
values (dots) and the theoretical expression in Eq. (|50p (con- 
tinuous line) of pf'^i (W)/p tl ^A(—W) for a 2D hard disk of 
mass Mo = 2 amu in a bath composed of 3 hard disks of the 
same mass. The initial energy is U = 0.831447 kJ/mol and 
the protocol doubles the mass of the disk. 



counts the number of occurrences of energy in the inter- 
vals J„ = [E n - AE/2, E n + AE/2) for a certain AE (in 
our simulations, typically, AE = 0.1 kJ/mol, for a total 
of about 20 intervals, and the histogram counts a total 
of about 10 5 events). This provides us with the start- 
ing statistics. At this point, we note that, independent 
of the functional form of M(t), acting the protocol on a 
particle with energy E gives with probability 1 the work 
W = E(M - M f )/Mf. The reason is that the time 
dependent system Hamiltonian {p\ +Py)/(2M(t)) gener- 
ates the following equation of motion for the momenta: 
p x = p y = 0. Hence E(t f ) = (jo? x + p 2 y )/(2M(t f )) = 
E(to)Mo/Mf regardless of the details of the protocol. 
So we immediately obtain a count of work belonging 
to the intervals J n = [W n - AW,W n + AW), where 
W n = E n (M -M f )/Mf and AW = AE(M -M f )/M f . 
After proper normalization, this yields a histogram, la- 
beled 



where we use the fact that the momentum (p x ,p y ) is a 
constant of motion. By applying the change of variable 
E = (pl+pl)/(2M ), and employing Eq. @J we obtain: 



l - 



M 



2 f 



W - U) /(CU) 



M — Mf 

Similarly one finds the backward pdf of work 

G,Uf, TI ,N rr-lr-, , ^-ll-<7 M 



C-l 



(48) 



Pt;7 J {-w) = uj l \i + c- 



1 - 



Ma 



M - M f 



\M f -M \ 

W - uA /(CU f ) 



c-i 



(49) 



Taking the ratio of Eq. 

p?; u t0 (w) 



and Eq. (f4"9"]) we obtain: 



By comparison with Eq. we see that the Finite Bath 
Fluctuation Theorem of Eq. ([28]) is satisfied. 



that provides a numerical estimate for Ptft W)- Next, 
for each n, we simulate the motion of the disk with fixed 
parameters, Mf = 2M and U n = U+ W n /{C + 1), and 
compute n different histograms for the backward prob- 



abilities hf g 'j™(k) in the same way as the forward his- 
togram was computed. By selecting the k — n value from 
each of the backward histograms and collecting them to 
form the new histogram 



c,u, 
h ,t 



we obtain a numerical estimate for Pt ' t f(— W). Finally, 

we compute the ratios h t ^ u to (n) /h to 'j™ (n) . 

These ratios are depicted in Figure [2] along with the 
theoretical values given by Eq. (f5CT)) . The figure shows 
excellent agreement between analytical theory and nu- 
merical experiment. The visible differences arc within 
the statistical errors. Note that, for the forward proto- 
col, where the mass is increased by a factor 2, the work 
can only be negative, and vice-versa for the backward 
protocol. Therefore, the graph shows only the negative 
values of noncquilibrium work W. 



V. DISCUSSION 



Numerical check of the Finite Bath Fluctuation 
Theorem 

In order to check numerically the validity of Eq. ([50]) 
we simulated the forward work pdipf- u to {W) for a bath of 
n = C 2D disks, a given value of U and a protocol that 
changes the mass of the disk from Mo to Mf = 2Mq. 
The pdf for the numerical work is calculated as follows. 
We first run a simulation of the motion of the disk with 
fixed U and Mq. We then construct a histogram that 



Physical meaning of 5T 

The basic quantity that enters the Finite Bath Fluctu- 
ation Theorem, and marks a distinction with the canon- 
ical fluctuation theorem of Crooks f|35[) . is the quantity 
ST, defined formally as the solution of Eq. (|2l)]) . This 
quantity enters in the definition of Uf and Tf . What 
is the physical meaning of these quantities? The hard 
sphere gas example turns useful in addressing this ques- 
tion. Calculations analogous to those leading to Eq. (14^]) 
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show that for a gas of hard spheres with a total of s de- 
grees of freedom, in contact with a bath with a specific 
heat C, it is: 



ST = W/C t , 



(51) 



where Ctot is the total specific heat of the system+bath 
compound system: Ctot '■= s/2 + C. This ST is therefore 
the increment of temperature that would result if, after 
having injected the energy W in the system of interest 
this is brought back into contact with the bath, and the 
compound system is let reach thermal equilibrium. Re- 
call that during the forcing protocol we assumed that sys- 
tem and bath are decoupled. We shall refer to this pro- 
cess as to the re-thermalization. After system and bath 
have re-thermalized, the extra energy W, initially stored 
in the system, will be shared between system and bath 
according to the ratio of the respective specific heats. In 
particular the bath gets the energy Q = CST, which is 
indeed the heat that flows from the system to the bath 
during re-thermalization. Accordingly the system looses 
this amount of energy, and its change in energy becomes 
AU = W — Q, in agreement with the first law of ther- 
modynamics. This means that U / represents the average 
energy of the system after the re-thermalization. To sum- 
marize: (a) The system is first in thermal contact with 
the bath. Its average energy is E/j and the temperature 
is Tj. (b) the system is decoupled from the bath and the 
forcing protocol is acted on it. As a result, the energy 
W is injected in the system with a certain probability 
density p t Jj (W). (c) The system (carrying the extra 
energy W), and bath (still at temperature Tj) are now 
allowed to rc-thcrmalize. During rc-thcrmalization the 
heat CST flows in the bath, the system reaches the av- 
erage energy U /, and the new temperature Tf is reached 
in the compound system. 

Remarkably, the temperature change ST vanishes in 
the canonical case: Hmc^oc ST = 0. However it is 
limc^oo CST = W , meaning that the whole extra en- 
ergy W injected in the system, flows into the bath during 
re-thermalization. However this does not affect its tem- 
perature (i.e., Tj = Tf), the specific heat being infinite in 
the canonical case. Therefore the term Tj/Tq does not 
appear in the canonical fluctuation theorem of Crooks. In 
fact the latter gives information about the free energy dif- 
ference of two states with different parameter values, but 
same temperature. This is a much more fortunate situ- 
ation as compared to the finite bath and microcanonical 
fluctuation theorems, in the sense that, in the canonical 
case, one should not bother to start the backward pro- 
cess from the "target" temperature Tf (which depends 
on W), but simply starts it from the same temperature 
as that of the forward process. 



From the canonical fluctuation theorem of Crooks, one 
obtains, after proper algebraic manipulations, and inte- 



gration over W, the integral form of the fluctuation theo- 
rem, namely the Jarzynski equality (e~P w ) = e~@ AF 0], 
which implies the second law in the form (W) > AF. A 
similar integral equation can be obtained for the Finite 
Bath Fluctuation Theorem too. It reads: 



N{T. 



•C-l e Sf{U f ) 



rpC-1 S {U ) 



where, 



M:= I p%%(W)dW 



(52) 



(53) 



Implications for the second law of thermodynamics rem 



and (•) denotes average over the normalized distribution 

q to , tf {W) := p U tlt s {W)/N. Eq. §5$ generalizes both 
the canonical Jarzynski equality and the microcanonical 
Entropy-from-work theorem [l3l . |27| |. Note that, as for 
the Entropy-from-work theorem, in general it is M ^ 1 
because the energy U / in Eq. ([53]) is a function of W (see 
Eq. As pointed out already in [23|, this prevents 

obtaining the second law directly from the integral form 
of the fluctuation theorem. 

Nevertheless the validity of the second law of thermo- 
dynamics for a driving protocol acting on a system that 
is initially thermalized with a finite bath, can be proved 
directly without invoking the Finite Bath Fluctuation 
Theorem. To this end it is sufficient to recall the content 
of two theorems which have been recently reported in 
the literature [H, [2{| [3(| ■ According to these theorems, 
the second law of thermodynamics, in either the minimal 
work principle form, or the entropy increase form of Clau- 
sius, is obeyed whenever the initial phase space pdf p(z) 
is a decreasing function of energy, namely p(z) > p(z'), 
for every z, z' such that H(z) < H(z'). This condition 
is obeyed by the Finite Bath Statistics, if the condition 
C > 1 is met (see Eq. In this regard we notice that 
this condition only is violated in the extremal case when 
the bath consists of a single degree of freedom (in which 
case it is C = 1/2), or if there is no bath at all (C = 0, 
microcanonical case) . 

The Crooks fluctuation theorem ([33)) can be seen as a 
statement according to which the probability of doing a 
certain negative work ~W during the backward protocol 
is exponentially suppressed with respect to the probabil- 
ity of doing the positive work W , in the forward protocol. 
For a cyclic protocol, this says that it is exponentially 
more probable to spend energy, rather harvesting it, in 
agreement with the Kelvin postulate (i.e., no energy ex- 
traction from a cyclic process). A similar situation oc- 
curs for the Finite Bath Fluctuation Theorem, with the 
exponential suppression being replaced by a power-law 
suppression. To exemplify this, consider again the gas of 
TV hard spheres in d dimensions. Imagine the protocol 
consists of changing the volume of the box that contains 
the gas from Vq to Vf. Straightforward calculations lead 
the following form of the Finite Bath Fluctuation Thco- 



p?; u t0 (w) 

p?:S(-w) 



vj 



w 



CtotTo 



c tot 



(54) 
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where it is evident that the power-law term (1 + 
W/ (CtotTo)) * ^ 1 becomes the exponential term ap- 
pearing in the Crooks theorem (|35|) for very large C 
{Ctot = C + dN/2 becomes very large for very large C). 



VI. CONCLUSIONS 

We devised a Finite Bath Fluctuation Theorem that 
gives information about the probability of work on sys- 
tems that have been thermalizcd with a finite heat bath. 
This corresponds to physical situations which are situ- 
ated between the two ideal cases of absent bath (micro- 
canonical ensemble) and infinite bath (canonical ensem- 
ble). The Finite Bath Fluctuation Theorem interpolates 
between microcanonical and canonical fluctuation theo- 
rems. It thus generalizes these theorems and reveals a 
common underlying mathematical structure. 

The validity of the Finite Bath Statistics is illustrated 
by means of numerical simulations of a 2D gas of hard 
disks in a box with perfectly reflecting walls, see Figure[TJ 
and the validity of the Finite Bath Fluctuation Theorem 
is confirmed both analytically and numerically, cf. Figure 
for our system. 

Similarity and differences between the Finite Bath 
Fluctuation Theorem and the canonical and microcanon- 
ical fluctuation theorems have been discussed, as well as 
its interrelation with the second law of thermodynamics. 
In contrast with the canonical fluctuation theorem, two 
temperatures, instead of one, appear in the Finite Bath 
Fluctuation Theorem. The physical meaning of these 
two temperatures has been clarified by considering a re- 
thermalization process. 

As shown in Sec. HU the Finite Bath Statistics in (|BJ) 
is a special instance of the general statistical formula ac- 
cording to which the bath density of states determines 
the shape of the system pdf. Based on quasi-adiabatic 



perturbation theory of chaotic systems, Jarzynski 31 1 
found that a slow particle coupled to a small bath with 
fast chaotic degrees of freedom thermalizes and reaches a 
stationary pdf whose shape is dictated by the density of 
states of the bath. Our simulations provide an example 
that such behavior of the system pdf occurs even if there 
is no time-scale separation between system and bath. In 
any case, thermalization of the subsystem towards a pdf 
of the form in Eq. ^ is expected only if the total system 
is ergodic. 

An important assumption underlying our main finding 
is that we used a specific heat that is energy- independent: 
Whether a Finite Bath Fluctuation Theorem exists also 
in the case of more realistic energy dependent specific 
heats remains an open challenge. 
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APPENDIX A: SPECIFIC HEAT OF A BATH OF 
n HARD SPHERES 

Although straightforward, the calculation of the mi- 
crocanonical specific heat of a gas of hard spheres is not 
discussed in statistical mechanics textbooks. We present 
this calculation below. 

The Hamiltonian of a gas of n d-dimensional hard 
spheres of radius a reads: 

n — >2 

H B ({pl}: {ql}) = E + E V (\Qi £1), ( A1 ) 



i=i 



where pi, ~ql are the d-dimensional momentum and posi- 
tion vectors of the i th sphere, and 



V(x) 



x > a 
+oo x < a 



(A2) 



is the hard core interaction potential. The phase space 
volume <&b with energy below Eg becomes: 

~ n „ n 

<s>b{e B ) = j / 

■—\ i—i 

n — >2 



i=i 



Kj 



(A3) 



where each integral in dqi is restricted to the region V, 
of volume V, of the box. For values of | ~ql — qj\ smaller 
than a, the integrand vanishes, thus reducing the spa- 
tial integration domain to the region M. C V™ where 
\qt — qj\ > a, for each couple In this region the 
interaction term is zero and one obtains: 



*b(E b ) =V 



f n ( n — >2 \ 

yn«»(*--E£:). 

i—1 \ 4—1 / 



(A4) 



where V' n = fuYl7=i^^i^ ^ s independent of Eg- We 
shall refer to V' as to the reduced volume. The integra- 
tion over the momenta then yields [32| : 

<5> B (E B ) = A dn {2m) dn / 2 V' n E% /2 (A5) 

where A N := n N / 2 /T(N/2 + 1). By differentiating 
^b{Eb) with respect to Eb, one finally obtains the den- 
sity of states of the gas of hard spheres: 



n B (E B ) = A dn {dn/2){2m) dn / 2 V lN E^' 2 - 1 



(A6) 



The only difference with the density of states of an 
ideal gas is that the actual volume V is replaced by 
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the reduced volume V . The temperature Tb(Eb) = 
$b{Eb)/Qb{Eb)i is given by same formula as for the 
ideal gas, i.e., Tb(Eb) = 2£ l s/(dn), and so is the spe- 
cific heat, i.e., C(Eb) = dn/2. For simplicity, in Eq. 
(|A1[) we neglected the spheres rotational degrees of free- 
dom. These however would add to the total specific heat 
an energy independent contribution. 



where we neglected the terms U as compared to CT. By 
making the change of variable x = CT—e, and neglecting 
again the term U as compared to CT, we obtain: 



h(U,T) 



/ dxn x (CT - x){CT - x)x°- 1 (B6) 
Jo 



APPENDIX B: EXISTENCE AND 
(NON)UNIQUENESS OF SOLUTIONS OF EQ. © 

We prove that, given U and A, it is always possible to 
find a T such that Eq. ([5]) is satisfied. For this purpousc 
we define the function: 



c-i 



pCT+U 

h{U,T):= deQ x (e)(e-U)(CT -e + U) 

Jo 

(Bl) 

which is continuous with respect to both U and T. The 
symbol il\ (e) denotes the density of states of the Hamil- 
tonian iJ(z,A). Eq. ([5|) can be equivalently expressed 
as: 



h(U,T) = 



(B2) 



For T = it is: 



h (U, 0) = / deSlx (e)(e-U)(U-e) 



c-i 



(B3) 

Since £l\(e) > 0, and e — U < in the integration domain, 
we have 

h(U,0)<0 (B4) 
On the other hand for T > U/C, we find 



CT 



I\(U,T)~ / den x {e){e-U){CT~er- 1 (B5) 
Jo 



All three terms forming the integrand are nonnegative, 
hence: 



h{U,T^>U/C) > 



(B7) 



Thus I\(U,T) is nonpositivc for T = and nonnegative 
for very large T . This implies, that there must be at 
least one nonnegative value of T, for which I\(U,T) = 0. 
Uniqueness, however is not guaranteed. 

In a similar way it is also possible to prove that 

h(0,T)>0 , I x (U^>CT,T)<0 (B8) 
showing that one can also fix T and find a U such that 
I\(U,T) = 0. Also in this case only existence is guaran- 
teed but not uniqueness. 

Examples for which two or more different energies 
correspond to the same temperature were reported in 
[33l . 134 ] for microcanonical (C = 0) gases with inter- 
particle interaction of the Lcnnard- Jones type. These 
systems undergo a microcanonical phase transition whose 
signature is the appearance of oscillations in the function 
T(U), which, therefore, is not invertible (i.e, U{T) is mul- 
tivalued). These oscillations are expected to appear also 
if these Lennard- Jones type systems are thermalized by 
means of a finite bath with specific heat C > 0. Based on 
the observation that no oscillation appear in the canon- 
ical treatment |34j |. one expects that the amplitude of 
these oscillations decreases with increasing C. 
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